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Formulas are obtained for computing the buckling load 
of rods eccentrically loaded at each end, the computation 
being extended in particular to the inelastic range. The 
test results are graphically presented on three sets of 
curves. Two of these, at least for the elastic range, are 
independent of the material tested. The third set (see 
charts 3a and 3b), which is independent of the material, 
possesses greater clearness and is therefore used for com- 
paring the test results with the thoorotical. On chart 6 
a comparison is made between the Duckling load of an ec- 
centrically loaded opon profile and the torsional buckling 
load of the same profile centrally loaded. Eor large slen- 
derness ratios the eccentrically loaded rod can sustain a 
greater load in the axial direction. 



NOTATION 



cross-sectional area of rod 

distance of outer fibers, with stress <J Q , 
from neutral axis 

a 

moment of inertia, J = i f 
radius of gyration* 
resistance moment 
eccentricity (fig. l) 
see. fig. 1 

slenderness ratio of rod 
deflection at center of rod 
deflection of neutral axis, v = f(u) 



f (cm 3 ) 
a (cm) 

J (cm 4 ) 
i (cm) 
W (cm 3 ) 
e (cm) 
u, l(cm) 
A = l/i 
8m (cm) 
v (cm) 



* "Knickbiegevorsuche. " Luf tf ahrtf or schung, vol. 17, no. 10, 
October 26, 1940, pp. 306-313. 
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E 


(kg/cm s ) 


elasticity modulus. 


T 


(kg/cm 2 ) 


buckling modulus, according to Engesser 


p 


(kg) 


axial load 


M 


(cm kg) 


moment, a function of the axial load P 


P Br 


(kg) 


axial load in failure 




(kg) 


bur*kl in.f InaA* P— — tt 3 -P f /\ . 




(kg) 


.2 

Euler load: P« = rr a fE/X 


E = 


. P E 
P 




X - 


e a 




7 - 


- ?Br _ °Br „ ^Br 


r* - 


. °0 _ *0 


,3 

— c - °° 


o — 


°¥ ~ tt 2 


~ » c o - 5= — 
I E 0 


cr 


(kg/cm 2 ) 


stress 


°0 


(kg/cm s ) 


stress in outer fibers in failure from 
bending and axial load 


°e 


(kg/cm 3 ) 


buckling stress: cr-g = ■n 3 l/X 2 




(kg/cm 2 ) 


buckling stress of Euler load: cr_ = tt 2 E/^ 

E o 


a 

g 


(kg/cm 2 ) 


limiting stress for which the Euler formula 
ceases to be valid 


a Br 


(kg/cm 2 ) 


axial stress at center of gravity during 
failure 


^e 


(kg/cm 2 ) 


bending stress (fig. 2) 
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■ 1. IKTR0DUCTI01T • 



There are several problems in connection with eccen- 
trically loaded beam columns, for which a theoretical solu- 
tion, is as yet unavailable. Although the problem may be 
theoretically formulated, it is necessary to .rely on tests 
for an explanation* and final solution* Two such problems 
aro the following: 

1. In eccentrically loadod boam columns, what impor- 

tance should bo assigned to the buckling mod- 
ulus ? 

2. lor which total stress P/f + M(P)/W is the beam 

no longer able to support any increased load? 

For the purpose of answering these questions the Focke- 
Wulf firm has undertaken a number of tests on beams loadod 
as shown in figure 1, 



,2. COMPUTATIONAL FORMULAS 



The beam deflection was- computed on the principlo 
of "passive work." The deflection at the center 5 m is 
obtained from the computation of the passive work 1 5 m 
performed by the assumed load "1" - this load acting at 
the location in the direction of the required deflection. 4 - 
when the actual load. P is applied to the beam. Prom the 
known relation, there is obtained: 



I /2 



i 6 „ = 2 f: -iL^£lJLUk^l du 



m 



E J 



1/2 



c -n ^ 2 / P v u . 

5 m = P o r~=rr + / du 



8 EJ . / E J" 
o 



The term is explained in greater detail in the author's 
paper "Tafel einiger Knick-Biege-Falle , " which will 
appear shortly. 
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For tho deflection v, we assume an oppression which 
satisfies tho boundary conditions for u = 0; u = 1/2 
and u = I : 

v = 6 m sin(TT S) 
There is then obtained: 



i 2 



a e j 



The maximum bending moment acts at the center and can 
be computed from the following equation: 

M max = ? e + p 5 m 

After several transformations, there is obtained for 
M max the following expression: 

t; to 1+ 0 . 233 , n x 

l£ max = p 0 g „ 1 — (1) 

The above is true on tho assumption that the buckling mod- 
ulus is equal to the S modulus. It will be explained be- 
low how this relation applies also in the inelastic range. 

We now have an expression for the maximum stress: 

cj _ I + 1.+ 0...233 < . v 

°max *" f + w |-1 - 0 1 ; 

where a can approximately be set equal to the compres- 
sive strength of the material. For thin-walled sections 
it is advisable to substitute for. this the "collapsing" 
strength - by which term we may denote the longitudinal 
stress obtained from compression tests on short-length 
specimens, of slenderness ratio < 10. 

An answer has thus been found to the second question 
in the Introduction. Tho theoretical answer, however, can 
only be confirmed by comparison with the test results. 

To express the results of equation (2) in the form of 
a curve chart, the relation is transformed for the limit- 
ing case of maximum load so that <y mSL x = <J 0 . There is ob- 
tained 
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'.' •- + 'y - .(i. + c) 

— . (3) 



1 + 0.230. y 
where 



oy "" P Br' °Br " e 



3. GRAPHICAL PRESENTATION OP THE RESULTS 



Chart- 1 shows the graph of equation (3). The abscissa 
is the 'aondimensional ratio x = e a/i s , which is an ex- 
pression for tine eccentricity. The following device was 
made use of in plotting. The curves we.re plotted for val- 
ues of x from 0 to 1. Por x > 1 , the reciprocal val- 
ues of x were taken and thus y as a function of l/x 
from 1 to 0. In this way all possible value's of x are 
included in the figure. The magnitude c 0 , which includes 

the maximum stress, the slenderness ratio, and the E modu- 
lus or buckling modulus is tatem as; the parameter of the 
curves. The ordinate is the quotient y Q = P^ r /Pj; o • This 

ratio is equivalent to p B r / p S' provided the longitudinal 
stress is small enough so that the computation with the 
elasticity modulus E is valid. In order to take into 
account also short rods of small eccentricity in which 
longitudinal stresses P/f may occur which exceed the 
stress o"gi the rupturing load P-g r for these rods was 

also referred to the Eulor load computed with the E modu- 
lus. This, in addition to computational advantages, has 
the following advantage for tho representation of the re- 
sults: If P Br wore referred to ' P-g, then all' the c 0 

parameter curves for x = 0 would end at point y = 1. 
This would mean large inaccuracies in the family of curves 
for small values of' ' x, and c Q •'• 

On chart 1 is also shown the limiting curve. To the 
right of this curve the longitudinal stresses are so small 
that the buckling computation with the E modulus is ab- 
solutely justified. To,, the. left of the limiting curve, 
however, the trend of the family of curves c Q depends on 
the material, the position of" the limiting curve itself 
depending on the material and to somo extent on the shape 
of the rod cross section. 



s 
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To determine the limiting curve, an assumption, is made 
which is fundamental as regards question 1 in the Introduc- 
tion. It is assumed that the longitudinal stress of the 
"beam P/f determines whether the "buckling computation, is 
to be carried out with the B modulus or with the buck- 
ling modulus. If P/f > o" , then the buckling modulus 

corresponding to tho state of stress is to ho used in the 
computation; while for P/f < Og the E modulus is to he 

used. The stress cr_ is thus the criterion in determin- 

§ 

ing whether the E modulus or the buckling modulus is to 
be used. 

The method is clarified in figure 2. Assume a rod is 
to be computed of slondernoss ratio X x . Dopcnding on the 

eccentricity, the rod can take up longitudinal strossos 
which arc smaller or greater than Og, up to tho failure 

load. If they are smaller than 0"_ the E modulus is to 

be used. For the slenderness ratio , the limiting 
point is 

Br g g 

_. = y g = c 

E 0 X E o x 0 

where 0" g and 0" 0 depend mainly on the material but also 
on the shape of the rod. In the relation 

(To. 

y g = t% o (4) 

we now have a simple expression for computing y as a 

function of c. Substituting (4) in equation (3) in order 
to make c approach zero, there is obtained: 

Ac 1 — y^ 
x = _ g ( 4a ) 

CT g 1 + 0.233 y 



or, solving for y 



Act 

Og - * 

s- + 0.233 x 
g 
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To compute the family of curves with parameter c 0 , 
a) an expression for the buckling modulus i s . introduced 
in equation (3), namely: 

cr 0 - a 
T = E — A 

Act 

and there is thus obtained 



Act 



^o/ 0 " - 1 ~ W. 



'O 

J o 

1 + 0.233 ~ G ° 



°o ct 0 /ct - 1 



Our object, however, is to express x as a function 
of ^Br/^Bo* Hence the expression takes the following 
form : 

, . Act 
c/y - 1 - -gr- c 0 

x = _2 , ( 5 ) 

1 + 0.233 ~ 

Relations (4a) and (4b) can, of course, "be derived 
from equation (5) with the aid of (4). It should also "be 
remarked that c/y = c 0 /y 0 . 

b) An expression for the buckling modulus according to the 
Tetmajer lines is introduced in equation (3), namely: 

" ^ W/Ag- J 

There is then obtained the following expression for 
the family of c D parameter curves to the loft of the lim- 
iting curve: 

Act 3 ct„ 

o- fa _ l £ 



(CT Q - CT) CT CT 

x — , ~_ « 

Act 2 o_ 
1 + 0.233 *S 
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Again transforming the abovo so that x is a functions 
sro is obta: 

(c/y - 1) - 



of 0"B r /°ii » there is obtained 



x = 




(6) 



In chart 2 the ordinate is the same as in chart 1, but 
for the abscissa there is taken the ratio o Q = °o/°]!! 0 

with x as parameter. The limiting curve can be very sim- 
ply drawn in by (4). To the right of the limiting curve 
the B modulus is to be used for the computation, while 
to the left of the curve the buckling modulus is used. 

Still a third form of graphically representing the 
equations (3), (5), and (S) will now be shown in chart 3. 
It is clear that these curves can be drawn using the val- 
ues for charts 1 and 2. The relations will, however, be 
modified to correspond to the new forms of representation. 
We are required to find cr-g r as a function of A with x 
as parameter. After several transformations there is ob- 
tained the relation: 

°o/°Br " f 1 " °' 2S3 *> 

In the above expression the significance of the assump- 
tion as regards the use of the buckling modulus for beam- 
column computations is particularly well brought out. The 
limiting curve in this case is simply the horizontal 0" = 
0" . Above this line relation (7) must be transformed and 

becomes : 

a) with the expression for the buckling modulus 




A 2 = 33 ™* Op/c^r I g o " °"3r (g) 

a 3r "o^Br - x) Act 

b) and with the buckling modulus according to Tetmajer, 
straight lines: 
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2 



T 




cr 



Br 



2 



a °./%r - (1-0-233 x) 



°o/°Br - ( 1+3C > 



v o' Br 



1) 




(9) 



By means of the family of curves with parameter x, 
charts 3a and 3b, the stress due to the axial loads on the 
rod is represented. The total stress which can be sus- 
tained by the rod in the outer fibers at the middle of its 
length is, according to the definition, cr 0 . Those rela- 
tions arc again presented in figure 3. For the cccontric- 
ity x the rod can take up an axial load corresponding 
to the stress °B r . The stress duo to the bending is 
then cr e . 

To show more clearly the conditions on increasing 
the axial load, the outer fiber stress cr = P/f + M(P)/W 
is plotted in figure 4 against the axial load P for var- 
ious eccentricities for a steel tube of diameter 30 X 1 
and of length, 100 centimeters. The dotted line gives the 
increase in purely axial stress with increasing axial load 
that is, for eccentricity zero. The buckling load P s = 
1890 kilograms. It may also be seen that even if greater 
stresses than the Bui or stresses in the outer fibers are 
permitted, no greater longitudinal forces than those of 
the Buler load can bo computed. Por the eccentricity of 5 
centimeters .there is also shown dotted the curve P/f + 
M/W without buckling deflection. Por this large eccen- 
tricity, therefore, the buckling effect is very small. 
It consists of the difference between the continuous and 
dotted curves. T7ith decreasing eccentricity the curves 
approach the linear stress p/f and the limiting curve P 

If the end points of the family of curves for the 
stress of 6000 kg/cm 2 , namely, the values for P Br were 
plotted against the eccentricity, a diagram similar to 
chart 1 would be obtained. The example here chosen corre- 
sponds approximately to a value of c = 3. 

4. COMPARISON OP THE THEORY WITE TEST RESULTS 



TTith regard to the employment of the buckling modulus 
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and the maximum stress, the procedure now is the following: 
Prom the tests, certain values x, c, A, etc., are deter- 
mined. In charts 4, 5, 6, only those x-parameter curves 
were plotted for which tests were conducted. Charts 4 to 6 
arc analogous to chart 3, sinco tho tost results are most 
clearly represented in this form. Comparison with the test 
values provides a critique of the validity of the proposed 
computation, methods. 

The results are typical for Duckling tests, particu- 
larly for small slenderness ratios. If such tests are not 
conducted with laboratory accuracy there is a large scat- 
tering of the results. The tests conducted cannot, however, 
be considered as laboratory tests and in practice, buckling 
members are not generally manufactured to laboratory pre- 
cision. With regard to each of the charts 4, 5, and 6, 
the following remarks may be made* 

Chart 4: The points with parameter value x = 2 lie 
quite well on the curve and similarly, tho points with 
parameter x = 1.33 and x = 1. Only for the values x = 
0.5 for slendernoss ratio A = 40 is there a strong scat- 
tering. Comparing, however, the values x = 0.5 for A = 
40 and A - - 70, it is found that there was hardly any 
increase in the failure axial load on decreasing the 
length from a value corresponding to the slenderness ratio 
70 to that corresponding to slenderness ratio 40. As this 
result is improbable, these points should be checked anew. 
Tho two upper points x = 0,5 and A = 40 indicate that 
the value of cr = 3400 kg/cm 2 has also been chosen some- 

O 

what high. The table below gives the ratio of the eccen- 
tricity to the mean radius of the circular tube used in 
the tests. 



x — > 


0.5 


1 


1.33 


2 


e / r m 


0.25 


0.5 


0.67 


1 



Chart 5: The test points lie above the theoretical 
curve. They could be made to lie more closely on the 
curves by taking cr somewhat larger than 3800. Tho 

curve x = 0 indicates the need for caution, however, in 
assuming a higher value for cr . 



Chart 8 : This chart shows tests with C sections 
connected as shown in figure 5, The tost points lie well 
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on the theoretical curves. The "buckling curve of a cen- 
trally loaded team which is capable of being freely twist- 
ed about its shear center is also shown for comparison - 
the curve being taken from the work of Eappus, Jahrbuch 
1937 der Deutschen Luf tf ahrtf orschung, page I 409. It will 
here be seen that for the usual slendornoss ratios the ec- 
centrically loaded section can sustain a greater longitudi- 
nal load than the centrally loaded section which has a less 
favorable degroo of froodom. Tho computation of the tor- 
sional buckling values is given in the appendix. 

The eccentric buckling tosts wero oxtonded by joining 
those and similar rods into a framework. In this case tho 
buckling value's of the joined rods are naturally somewhat 
higher. This is principally due to their being joined to 
the other rods under tensile load. 



5. C01T0LUDIHG REMARKS 



On the basis of the tost values and from comparison 
of these values with the theoretical curves, it may be 
stated that the buckling process may, with sufficient ac- 
curacy, be represented by the theory given above. Exami- 
nation, of the curves leads to the following observations. 

It is preferable to conduct tests on eccentrically 
loaded rather than on centrally loaded beans. Generally 
great difficulty will be found in conducting buckling 
tests on centrally loaded beams. The difficulty consists 
in the fact that the buckling value of a centrally loaded 
beam is dot erminod by the offset of the neutral axis at 
each section from the line joining the points of applica- 
tion of tho load. Tho amount of this tolerance, which is 
not capable of being measured, is decreased by the pres- 
ence of an eccentricity. In buckling tests on centrally 
loaded beams failure occurs due to the moment P A where 
A is the tolerance. In compression tests on eccentrical- 
ly loaded beams, failure is due to the moment P (e + A). 
Erom this it may be seen that the effect of the tolerance 
vanishes with increasing value of e. 

The tests reported in the present' paper are only those 
conducted by Eocke-TCulf . Results of other provious tests 
were also evaluated (Thielemann, Jahrbuch 1937 der Deut- 
schen Luf tfahrtf orschung, p. I 386). The results agree 
almost completely with our own. One phenomenon, however, 
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is there brought out more clearly than in our paper, Tor 
large slenderness ratio the test points lie somewhat below 
the theoretical curve, while for smaller slenderness ratio 
they lie above - and the more so the greater the deflec- 
tion due to bonding (fig. 6); that is, tho greater c e . In 

other words, for large ' slenderness ratio the stability part 
predominates in tho failure of tho beam, and for small 
slenderness ratio the stress part predominates and an ef- 
fect arises that is called tho "bending factor." (See ref- 
erence 1.) 

Besides the references already given there may be men- 
tioned, from the very extensive literature on the subject, 
the book of Karl Jezek (reference 2) . This book gives a— 
very detailed treatment of the flow processes or what has 
boon noted above as "bending factor" and the effect of the 
cross-sectional shape. The book is intended particularly 
for steel structural problems. 



APPENDIX 

Torsional Buckling of Centrally Compressed Open Sections 
in the'Elasthc Eange (Reference 3) 



According to equations (47), (48), (5S) , and (58), we 
have for the torsional buckling stress cr^ of a section, 
symmetric with respect to the x-axis: 



(ii) i 



.°S + °3 




J 



\2 

S ^ + p 



TT 2 E J , 



I* 1 



TT' 



0"„ = 



& J T + p E C 



rr e a. 



J 1 J, 



where 
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r J x moment of inertia with respect to the x-axis 
polar moment of inertia 

■I g 

Jrp * U s , torsional resistance (U = perimeter, s = 

thickness) 

'S Ej. / y w & J, - deformation moment for tho x-axis 

0 / w d P - j (wdP), resistance to deformation about 

the center of gravity 

v w unit deformation 

Por a U-soction, equations (I) and (II), "by approxi- 
mating tho roots of equation (I), can Do put in tho form 
(see equations (69) to (72) in the above paper (reference 
3)): 



(IVH 



bj_| j*- (HI) cr D = Sf| |;a + p A 2 ] 

A i + As //A, - A a f " 

2 " / V 2 / + 



a = 



fi = 



3 
3 



1 (> " ^) 




'A x = 3a 3 (a + 2b) (2 a S + 15 ah + 26 b 3 ) 

(2a + b) [a 2 (a + 6b) (a + 2b) + 4b 2 (2a + b)] 



00 J 



A a = 



A* = 



B 



2 

a (a + 6b) (a + 2b) 
4b 3 (2a + b) 

6a (a+2b) (a-f3b) 

b ( 2a+b ) Jz[a? (a+6b ) (a+2b ) + 4b s (2a+b J] 

(a+2b) 2 s s 

[a S (a+6b) (a+2b) + 4b 3 (2a+b)] 
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In tho special cases there is obtained: 



20 X 20 X 0.5 

I I 

20 X 15 X 0.5 



TT 2 B 
A U 



tt" e r 

A L 



O.XB + 0.079 (±.) j 
(£o)] 



0.31 + 0.136 



A = 
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Figure 1.- Eccentrically 
loaded rod. 
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Figure 2.- Buckling stress 
diagram. 
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Figure 5.- Connection 
for test 
C section. 
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For comparison, torsional "buckling;. 
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holds for 
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67 



(l-0.233x) 6Br =<5 g 
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J curve 
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J curve 



, 6p-CBr holds for 
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E = 7 •10 5 kg/cm! 
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Chart 6. 
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Chart 1. 




y= £l 



0.4 0.8 0.8 0.4 
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"££±0.233x 
V 
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1+0.233 y 
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limiting curve & ~ <T " C ° 



valid to right 
of 

curve 



A6~. Cn 
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Chart 1. 
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Chart 2. 
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Chart 2. 
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